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Abstract 

In certain neutron-rich Te isotopes, a decrease in the energy of the first excited 2"*" state is accom- 
panied by a decrease in the E2 strength to that state from the ground state, contradicting simple 
systematics and general intuition about quadrupole collectivity. We use a separable quadrupole- 
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• plus-pairing Hamiltonian and the quasiparticle random phase approximation to calculate energies, 

B{E2,0~^ — > 2"*") strengths, and g factors for the lowest 2"*" states near ^'^^Sn {Z > 50). We trace 
the anomalous behavior in the Te isotopes to a reduced neutron pairing above the N = 82 magic 
gap. 
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I. INTRODUCTION 



As experiment pushes towards the nuclear drip hne, it is becoming possible to examine 
isotopic chains over increasingly large ranges of and Z . We have new opportunities 
to test systematics and the ideas that underlie them. One region in which experimental 
progress has been made recently surrounds the neutron-rich doubly magic isotope ^^^Sn. In 
particular, recent Ref. [0 reports measurements of the transition strengths B{E2] 0+ 2+) 
(or B{E2)^ for short) from the ground state to the lowest 2"*" state for ^'^^Te, ^^^Te, and 
^^^Te. The authors discovered that B{E2)ys and the energies of the lowest 2"*" states (-E2+) 
behave differently in the Te isotopes (with = 80, 82, and 84) than in those of Xe, Ba, 
and Ce which have more protons. In most isotopic chains, including those three, a decrease 
in E2+ is accompanied by an increase in B{E2)^ as the states become collective. This is not 
the case in i32,i36'jg^ where the B{E2)^ decreases as E2+ decreases. 

Our work explains this unusual behavior. Our tool is the quasiparticle random phase ap- 
proximation (QRPA), in conjunction with a simple schematic interaction, which we apply to 
even-even nuclei in the mass region 50 < Z < 58 and 80 < A^ < 84 (and a much larger range 
of A^ for the Sn chain). The QRPA is a well-established method for describing vibrational 
states PI and has advantages of simplicity, particularly when separable interactions are used 
and exchange terms neglected. One should mention, that there exist large-scale shell-model 
calculations for selected nuclei around ^'^^Sn [|I], ||, However, at the present stage, these 
calculations use different spaces (and interactions) for nuclei above and below the N=82 
magic gap. Our model, albeit more phenomenological, uses the same Hamiltonian in both 
regions. 

This paper is organized as follows: In Sec. II we review phenomenological and simple 
microscopic approaches to the systematics of E2+ and B{E2) |. In Sec. Ill we give an 
overview of the experimental data around ^^^Sn and discuss their significant properties. In 
Sec. IV we use the Hartree-Fock-Bogoliubov (HFB) method to discuss static properties of 
the ground states. The QRPA model is described in Sec. V. We show results of the QRPA 
calculation for the lowest 2~^ states in Sec. VI and discuss the origin of the irregular behavior 
of B{E2)^ from a microscopic point of view in Sec. VII. The g factors for Xe, Te, and Sn 
isotopes are treated in Sec. VIII. Finally, Sec. IX summarizes this work. 
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II. RELATION BETWEEN AND B(E2)'l 

The systematic relation between £'2+'s and -B(i?2)|'s is an old topic. One early phe- 
nomenological relation (by Grodzins 0) is 

B(E2;0^^2^)^14,9^-i^^[eVl, (1) 

and another (by Raman et al. 0) is 

B(E2;0+-.2-) = 3.26p-i^^leVl. (2) 

The latter reproduces most of more than 300 experimental data points to within a factor 
of 2. Both these formulae, after factoring out a gentle dependence on Z and A, assert that 
-B(i?2)|''s are inversely proportional to -E'2+'s- For vibrational states, this result is predicted, 
if mass parameter is constant, by the liquid drop model 0, which gives 

B(E2; = ^ = 1) = 5 {^^ZeR^)' (3) 

where R is the nuclear radius, and D2 the quadrupole mass parameter. 77-2 denotes the 
number of 2"*" phonons. It also falls out of an RPA treatment of collective excitations in 
the simple microscopic model of Brown and Bolsterli and others p|, 0. In physical 
terms, collectivity lowers the energy of attractive modes while at the same time increasing 
the transition strength because nucleons contribute coherently to the transition. 



Another successful way of classifying collective 2+ states is the NpNn scheme [11, |T^, 13 
Both the -E2+'s and -B(-E'2) fs lie on smooth curves when plotted as functions of NpNn, 
where Np (A„) is the number of valence proton (neutron) particles or holes. The plot for 
some nuclei around those considered in this work is shown in Fig. 1. The data points can 
be divided, approximately, into two well-correlated groups: those for < 82 (the upper 
E2+ and the lower B{E2)^ branches) and those for > 82 (the lower E2+ and the upper 
B{E2)^ branches). The plots reveal a clear asymmetry in the 2+ states with respect to 

= 82. That is, the > 82 systems have lower E2+ and higher B{E2)^ as compared to 
their A^ < 82 NpNn partners. This would suggest increased quadrupole collectivity in the 
region above A^ > 82. However, as discussed in the following, deviation from this general 
trend can be found. 
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FIG. 1: Lowest 2+ energies (top) and B{E2)]^s (bottom) versus NpNn in a number of even-even 
nuclei with 52 < Z < 64. The data are from Refs. [||, The curves are to guide the eye. 

III. OVERVIEW OF DATA AROUND ^'-^^SN (Z > 50) 

Let us survey the experimental data relevant to this paper. Figure 2 shows E2+^s and 
B{E2) I's for the lowest 2+ states of even-even nuclei as functions of neutron number. 
Both observables are fairly symmetric around = 82 for the Xe-Ce isotopes indicating 
that particle and hole excitations in those nuclei play similar roles. Actually, some of the 
B[E2)ys in Ce and Ba in the region > 82 are slightly larger than those with the same 
Nn in N < 82; similarly the E2+'s for > 82 are lower than those for < 82, in a way 
consistent with the NpNn plots of Fig. 1. Clearly these isotopes follow the usual relation 
between B{E2)^ and -£^2+- 

On the other hand, ^^^Te, ^^^Te, and ^^^Te behave differently. The B{E2)^ is not symmet- 
ric adjacent to A^ = 82, a fact that is even more significant when looking at the corresponding 
energies in Fig. 2. The state in ^^^Te lies 370 keV lower than that of ^^^Te, but nevertheless 
the B[E2) I in ^^^Te is smaller than that in ^^^Te. The situation violates the pattern of 
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FIG. 2: Experimental values of E2+ (top) and B{E2)'\ (bottom) in even-even Sn, Te, Xe, Ba, and 
Ce isotopes, as functions of neutron number A^. The experimental B{E2)] rates were taken from 
Refs. (for i?2+ cf. §). 



typical collective behavior discussed above. [This behavior does not appear anomalous on 
the NpNn plots of Fig. 1 because of the scale of the figure, however, NpNn = 4 for both ^^^Te 
and i36Te, and E2+(^^^Te) = 0.974 MeV, E2+(^^^Te) = 0.606 MeV, B{E2,^^^Te) = 0.172 
6%"^, and 5(E2,i36Te) = 0.103 e%\] 



IV. HFB CALCULATIONS 

As a prelude to our QRPA treatment of the 2+ vibrations, we calculate static shape and 
pairing deformations in the Hartree-Fock-Bogoliubov (HFB) model of Refs. fl^, [15|, |1^, [1^, 
p!8| . We perform axially deformed HFB calculations with the particle- hole Skyrme forces 
SLy4 ||19| and an intermediate contact delta pairing force ^7^. The resulting quadrupole 
deformation parameter /3 = J\\^Q) Q being total quadrupole moment and R - rms 
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FIG. 3: The quadrupole deformation parameter f3 calculated in the HFB approximation with the 
Skyrme force SLy4 and an intermediate-type delta pairing force [ p^ . 

radius, is shown in Fig. 3. It indicates that the static deformation of the nuclei with = 80 
and 84 is zero or small (~ 0.1) compared to those of the mid-shell nuclei. We can therefore 
treat the 2+ states in these nuclei as vibrations around a spherical shape. 

In general, the HFB calculations follow the NpNn trend discussed earlier. The (3 values 
above the = 82 gap are systematically increased for > 4. The strongest asymmetry 
in the pattern of (3 is predicted for the Te isotopes. 

Figure 4 shows predicted neutron pairing gaps. Since pairing is a symmetry-restoring 
interaction, the calculated pairing gaps are anticorrelated with the quadrupole deformations. 
Consequently, the values of A„ are systematically lower as one crosses the = 82 gap. In 
particular, in most cases A„(A^ = 80) > A„(A^ = 84). 
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FIG. 4: Same as Fig. 3 but for the neutron pairing gaps. 



V. QRPA CALCULATION 



The Hamiltonian we use in our QRPA calculation is 



(4) 



where is the single-particle energy, and cj^ is the creation operator of a nucleon in the 
state 11. Xr is the chemical potential, which depends on the isospin z-component r. A,- is 
the pairing gap, and is the monopole pair creation operator. 

As a residual two-body interaction, we use the sum of an isoscalar quadrupole force Hq, 
an isovector quadrupole force Hq, and a quadrupole pairing force Hq, defined as follows: 
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FIG. 5: The experimental s.p. spectrum of Sn (from |2. 
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where /i denotes the time-reversed state of /x. For Xt=0; we use the self-consistent values 
of Ref. for xt=i, we use the value xt=i = XT=i(std) = -92.9^-^/3 ^eV frn'^ (As 
will be seen later, the results of QRPA calculations are fairly insensitive to the choice of 
Xt=i-) We fix the quadrupole pairing strengths according to the prescription proposed 
in Ref. [We refer to this value as (self).] Our QRPA equations are in the standard 



matrix form, as in Ch. 14 of Ref. [22|, and, as usual, we neglect the exchange terms of the 
multipole-multipole interactions. 

Our calculations are performed in a single-particle (s.p.) space of several harmonic- 
oscillator shells {Nose = 2-6 for protons and A^^osc = 2-7 for neutrons). Since our configuration 
space is large, we use the bare, rather than effective, charges in calculating B{E2)^. We take 
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FIG. 6: The experimental neutron pairing gaps (connected by lines) obtained from the odd-even 
mass differences and calculated pairing gaps with the HFB-Lipkin-Nogami method (isolated sym- 



bols). Experimental masses are from Ref. [29| 



s.p. energies e:^ from experimental data around ^^^Sn, shown in Fig. 5. (When the levels 



are not available this way, we use Woods-Saxon energies p3[ for bound levels and Nilsson 
energies for unbound levels.) It is worth noting that the neutron level density just below 
the 82 shell gap is much larger than it is above the gap. This is due to the near-degeneracy 
of lhu/2, 2(^3/2, and 3si/2 shells and a fairly large energy gap between the 2/7/2 and 8^3/2 
shells. As we will see, this difference plays a crucial role in the anomalous behavior of the 
Te isotopes. 
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Figure 6 shows the experimental pairing gaps obtained from odd-even mass differences, 
according to the prescription of Ref. |2^, and gaps calculated by the HFB-Lipkin-Nogami 
method |2^. We note that the HFB-Lipkin-Nogami calculation, which partly corrects for 



particle number fluctuations, reproduces experimental trends very well. The neutron pairing 
gap in the Sn, Te, and Xe isotopes decreases as increases and crosses = 82. This effect, 
clearly seen also in the HFB calculation of Fig. 4, has been noticed earlier, cf. Ref. [|T7[. In our 



QRPA calculations, we used renormalized experimental pairing gaps. The renormalization 
factors, reflecting the reduction of pairing in excited 2^ states, were adjusted to experimental 
data in the Sn isotopes. The renormalization factor turned out to be 0.6 (0.9) for neutrons 
(protons). For magic nuclei with = 82 and/or Z = 50, we took A = 0.4 MeV, a somewhat 
arbitrary value, reflecting the weak pairing correlations in magic nuclei. (Experimental odd- 
even mass differences for magic nuclei do not determine pairing gaps well ||25| . ) We used the 
average of the proton pairing gaps at A^ = 80 and 84 for Ap at N = 82 to avoid the sudden 
decrease at the magic number. 

VI. RESULTS OF QRPA CALCULATIONS 



We carried out QRPA calculations for even-A^ isotopes of Sn with A^ = 64-84, and for 
the A^ = 80, 82, 84 isotopes of Te, Xe, Ba, and Ce, which are nearly spherical in our HFB 
calculations. Figure 7 shows the calculated lowest 2"*" energies and i?(£'2)|'s, along with 
the experimental data. The calculations reproduce the experimental trend quite well, in 
particular the asymmetry around A^ = 82 of the B{E2) |'s in the Te isotopes. We also 
predict an inverted, and more symmetric, curve for the B{E2)^^s in the Sn isotopes with A^ 
= 80-84. This kind of inversion is well known to occur in the Pb region around A^ = 126 [Q. 
(For more discussion of this point, see Sec. VII.) For comparison. Fig. 8 shows the results 
with the pure Nilsson spectrum (parameters from Ref. [Q). The collectivity in the A^ = 
68-76 isotopes of Sn is enhanced here, but otherwise Figs. 7 and 8 are fairly similar. Kubo 



et al. 1^ performed calculations in Sn isotopes up to A^ = 74 with a similar Hamiltonian 
and obtained similar results. In the shell-model calculation of Ref. [|l|, B{E2)\ for ^^^Te 
^i36rpg ^ turned out to be 0.088 (0.25) e^6^, i.e., the transition rate has been predicted to 
increase when going from A^=82 to A^=84. 
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FIG. 7: -E2+'s (left) and i?(£^2)|'s (right) from the QRPA calculation and the experimental data. 
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FIG. 10: Same as Fig. 9 but as functions of Gl/G^{se\i). 



We checked the stabiUty of our calculations by varying the strengths of the isovector 
quadrupole force and the quadrupole pairing force. Figures 9 and 10 show the results in 
Te. The unusual behavior around N — 82 clearly is not sensitive to the strengths of these 
forces. Based on all these results, we conclude that the QRPA prediction of the unusual 
behavior around ^^^Te is robust and does not depend significantly on model details, except 
for neutron pairing. 



VII. ABNORMAL PATTERN OF QUADRUPOLE COLLECTIVITY IN THE 
NEUTRON-RICH TE ISOTOPES 

What is the reason for the unusual behavior of the Te isotopes around N — 82, i.e. the 
fact that both E2+ and the B{E2)'[ are smaller in ^^^Te than in ^^^Xe? The ingredient in 
our calculations that displays the most asymmetry around A?" = 82 is the neutron pairing 
gap. To understand how it affects the results, we performed QRPA calculations in ^^^Te for 
different values of A„. The results are shown in Fig. 11. As A„ decreases from 0.6 MeV 
to 0.4 MeV, both the E2+ and B{E2)^ decrease, indeed suggesting that this quantity plays 
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FIG. 11: The lowest 2+ energy (top), B{E2) t (middle), and the summed QRPA amplitudes 

J^^vii^fiu ~ ^fi.u) protons and neutrons (bottom) as functions of the neutron pairing gap in 
i36'jig_ r^Yie arrows show the locations of the gaps in I32,i36'pg used in the solution in Fig. 7. 



the key role in the unusual trend we want to explain. To get more insight, we consider the 
forward {ipfj,u) and backward {(fi^u) QRPA amplitudes in the lowest-energy 2+ excitation 

|2^> = ^{i^i^^alal - (/?^^a^a^)|g.s.), (6) 

where aj^ and create and annihilate a quasiparticle in the state /i, and |g.s.) is the QRPA 
ground state. The QRPA amplitudes ipn^ and (pn^ depend on the ratios 
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i32Te i34Te ise^g i34xe i36xe ^^Sxe 



proton 
neutron 

E 'fill' 

proton 

E 'fil'^ 

neutron 



0.63 0.99 

0.44 0.02 

0.03 0.00 

0.04 0.01 



0.12 0.76 

0.97 0.40 

0.04 0.08 

0.05 0.08 



0.99 0.52 

0.04 0.67 

0.01 0.09 

0.02 0.10 



TABLE I: Summed squared forward (V'^,^) ^'^d backward (<^^,^) QRPA amplitudes for N = 80, 82, 
and 84 Te and Xe isotopes. 

respectively, where S^j, — ^Jis^ — A,-)^ + is the BCS quasiparticle energy. The bottom 
portion of Fig. 11 shows that these quantities depend significantly on the neutron pairing 
gap as well. 

The reason for the unusual behavior can be surmised from these figures. The decreased 
neutron pairing gap in ^^^Te means that the lowest neutron quasiparticle energies are lower 
than those in ^^^Xe (0.792 MeV for ^^^Xe and 0.460 MeV for ^^^Te). As a result, the energy 
of the lowest 2+ state decreases when one crosses N — S2. But the low- lying neutron 
quasiparticle energies also cause the neutron amplitudes in the wave function to increase 
and the proton amplitudes to decrease, as Fig. 11 and Tabic I show. Since the B{E2)] is 
determined solely by protons, it decreases as well. In other words, the behavior of the lowest 
2"*" states reflects properties of the s.p. spectrum — the fact that it is more dense below 
N — 82 than above (see Sec. V), giving rise to a larger pairing gap — more than collective 
quadrupole effects induced by the residual interaction. This is not a total surprise given 
that both isotopes have only 2 valence neutrons (or neutron holes). 

In the Xe, Ba, and Cc, isotopes, the increased number of protons makes proton pairing 
and the neutron-proton quadrupole-quadrupole interaction more important and reduces the 
effectiveness of the s.p. mechanism just described. This is nicely illustrated in Fig. 12 for 
^^^Xe. One can see the usual relation between E2+ and B{E2) f and a clear difference 
between Te and Xe in the A„-dependence of B{E2)^. In Xe, B{E2)^ increases as the proton 
amplitude decreases, indicating increased collectivity. 
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FIG. 12: Same as Fig. 11 but for ^^^Xe. The values of A„ in ^^^'^^^Xe, employed in QRPA 

calculations, are marked by arrows. 



The value of B{E2)^ in ^^^Te is smaller than that of ^^^Te, in spite of the large proton 
amphtude (see Table I). However, the 2+ state in ^^^Te corresponds to one two-quasiparticle 
configuration ((77/2)^, while the strength in ^^^Te and ^^^Te is more fragmented, indicating 
the collective character of the state. 

We close this section by discussing the behavior of B{E2)^ of ^^'^Sn-^^^Sn mentioned in 
Sect. VI (see Fig. 7). For this purpose. Fig. 13 shows summed QRPA amplitudes for protons 
and neutrons in the Sn isotopes. It is clear that the neutron amplitudes are dominant in 
all cases. However, at ^^^Sn, both proton and neutron low-energy excitations are hindered; 
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FIG. 13: Summed squared amplitudes J2ij,ui'^fiu ~'Pfiu) the protons and neutrons of Sn isotopes. 

therefore the neutron amphtude decreases and the proton contribution increases, compared 
to the other isotopes. This change causes a local increase in B{E2) | at ^^^Sn. (When 
the collectivity is small, B{E2)^ reflects the magnitude of the proton amplitudes directly.) 
Since the nucleus is in a neutron-rich region, however, matrix elements of the quadrupole 
operators of the neutrons are larger, on average, near the Fermi surface than those of the 
protons. Thus, excitations of the neutrons are still dominant in the 2"*" state of ^^^Sn. 



VIII. g FACTORS OF XE, TE, AND SN ISOTOPES 

The abnormal behavior of the -E2+'s and B{E2)^^s around ^^^Sn reflects the variations of 
proton and neutron amplitudes in the wave function of the lowest 2+ state. Therefore, we 
analyze the g factor in neighboring nuclei; they are very sensitive to relative proton/neutron 
compositions. 

We have calculated the g factors of ^^^Xe, ^^^Xe, and ^^^Xe, and compare with recent data 



| 27| in Table II. As usual, we multiplied the bare spin g^ factors by 0.7, and took bare gi 
factors |lOl. Our g factor in ^^^Xe is larger than in ^^^Xe, though not by as much as the 
data (see also Ref. 0). We show the corresponding proton and neutron QRPA amplitudes 
of 2+ states in Table I. Protons are more important in ^^^Xe and ^^^Xe, while neutrons are 
more important in ^'^^Xe. We found by analyzing the amplitudes that the main component of 
the 2+ states of ^^^Xe and ^^^Xe is 7^(1^7/2)^ while those of ^^sXe are 7i{lgy2y and 2/(2/7/2)2. 
It is interesting to compare the g factors with those of the single-particle states in Table III. 
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i34xe i36xe i38xe 



exp. 0.354(7) 0.766(45) 

cal. 0.585 0.716 0.291 



TABLE II: Experimental and calculated g factors for I34,i36,i38x;g isotopes. The data are from 
Ref. p7|. 



neutron proton 





fit 


th. 


fit 


th. 


2^3/2 


0.554 


0.534 


0.544 


0.419 


1^11/2 


-0.223 


-0.243 


1.39 


1.264 


3Sl/2 


-2.65 


-2.674 


4.04 


3.906 


2^5/2 


-0.514 


-0.535 


1.54 


1.581 


I57/2 


0.317 


0.297 


0.803 


0.677 



TABLE III: The g factors for neutron holes in ^^^Sn and proton particles in ^^^Sb. The values 
labeled as "fit" are taken from Ref. |27|, while the theoretical estimates are Schmidt values with 
gs multiplied by 0.7. 



132Tg 134xe 136xe 



0.491 0.695 -0.174 



TABLE IV: The calculated g factors of I32,i34,i36rpg iggtopes. 



The observed g factors for ^^^Xe and ^^^Xe support the idea that the states of these nuclei 
consist mainly of proton excitations (see Ref. |2^); our calculation is consistent with this 
picture. The large g factors of the proton 1/111/2, 3si/2, and 2^5/2 orbitals suggest that the 
nuclear g factors are sensitive to the small admixtures of these orbitals. The Xe isotopes 
therefore provide a severe test case of the many-body wave function. 

Table IV displays calculated g factors of the neutron-rich Te isotopes. The neutron 
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FIG. 14: Calculated (asterisks) and experimental [^0[ (open squares with error bars) g factors of 
the lowest 2+ states for Sn isotopes. 

dominance in our ^^^Te wave function clearly lowers the predicted g factor there. It would 
be interesting to test this prediction experimentally. 

Finally, Fig. 14 shows calculated g factors of Sn isotopes compared to the experimental 
data. The behavior of the g factors up to = 74 can be understood in terms of the negative 
single-neutron (yf-factors of the l/iii/2, 2(^5/2, and 3si/2 shells (see Table III and Ref. [pO[| ). 



Around N=78, however, the 2^3/2 orbital carrying a positive g factor becomes occupied, 
and this gives rise to positive ^f-factors in 128,130,132g^^ 

Above A^=82, the structure of the 
lowest 2+ state is dominated by 2/7/2 shell, and (/-factors drop again. 



IX. SUMMARY 



In this paper, we have investigated the irregular behavior of -E'2+'s and S(£'2) fs in 
^^^Te-^^^Te through the QRPA with a simple separable interaction. Our QRPA calcula- 
tions reproduce the behavior seen in experiment, and we trace the cause to the difference 
in neutron pairing below and above A^ = 82. The decrease in A„ with A^ is clearly seen in 
experimental systematics and in self-consistent calculations. The results of our phenomeno- 
logical model are fairly robust and depend only weakly on other model parameters. A related 
finding is that the B{E2)^ in ^'^^Sn should be larger than in its immediate Sn neighbors, as 
is the case around ^"^^Pb. We hope that this prediction will stimulate further measurements 
in the neutron-rich region around ^^^Sn. 
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To strengthen our argument about neutron dominance in the wave function of the 2+ 
state in ^^^Te, we also calculated g factors of the Xe, Te, and Sn isotopes. We reproduced 
the experimental trends and found that while protons dominate the excitation amplitudes 

in ^^^Xe and ^^^Xe, the g factor of the 2+ state of ^^^Te is dramatically reduced. The 
experimental discovery of this effect as well as significant behavior of ^^^Sn-^^"^Sn would 
validate our understanding of the structure of nuclei around ^^^Sn. 
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